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Chapter  I 
INTRODUCTION 

The  drastic  reduction  of  wall  shear  stress  induced  by  the  addition 
of  small  amounts  of  high  molecular  weight  polymer  to  a  solvent  flow  was 
first  observed  by  Toms  (1948)  twenty-four  years  ago. 

Since  then  many  theories  have  been  proposed  to  explain  the  drag 
reduction  phenomenon.  Toms  (1948)  postulated  the  existence  of  a  shear 
thinning  wall  layer  of  low  viscosity,  in  the  polymer  flow,  which  would 
yield  lower  friction  coefficients  than  in  the  solvent  flow.  However, 
rheograms  of  drag  reducing  polymer  solutions  obtained  by  conventional 
viscometry  methods  have  shown  Newtonian  rather  than  shear  thinning 
behavior. 

It  has  also  been  speculated  that  the  polymer  induces  an  anisotropic 
viscosity,  that  is,  a  viscosity  which  is  low  in  the  direction  of  flow  but 
high  enough  in  the  other  directions  to  dampen  turbulent  fluctuations. 
Differences  in  normal  stresses  may  provide  such  a  mechanism.  Gadd 
(1968)  measured  normal  stresses  in  several  polymer  solutions  and  found, 
with  one  exception,  that  there  was  no  significant  difference  in  the  stress 
components.  He  concluded,  therefore,  that  there  is  no  correlation 
between  difference  in  normal  stresses  and  effectiveness  in  drag  reduction. 

Lumley  (1967)  suggested  that  molecule  entanglements  extending 
across  the  stream  resist  the  formation  of  the  streamwise  vortices  and 
therefore  enlarge  the  sublayer  eddies.  In  a  similar  vein,  Peterlin  (1970) 
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sees  one  end  of  the  polymer  macromolecule  lying  in  the  core  of  a 
"microvortex",  the  other  end  lying  outside.  The  macromolecule  is 
thus  deformed;  while  being  deformed  it  absorbs  energy  and  dissipates 
the  vortex. 

Walsh  (1967)  has  proposed  one  of  the  more  interesting  and  compre¬ 
hensive  theories  to  date.  He  begins  by  assuming  that  the  disturbances 
responsible  for  Reynolds'  stresses  are  born  at  the  edge  of  the  sublayer 
and  that  the  primary  role  of  the  random  turbulence  in  the  wall  shear 
layer  is  to  excite  a  violent  repetitive  instability  in  the  sublayer  flow. 

The  breakdown  of  this  instability  produces  a  small  disturbance  at  the 
edge  of  the  viscous  sublayer.  This  small  disturbance  tends  to  grow 
because  it  can  locally  extract  energy  from  the  mean  velocity  profile 
through  momentum  transfer.  Simultaneously  this  small  disturbance 
tends  to  lose  energy  because  its  gradients  locally  dissipate  energy  into 
heat  and  locally  store  energy  in  the  polymer  molecules.  If  a  small 
disturbance  extracts  more  energy  locally  than  it  loses,  it  will  grow. 

The  disturbances,  or  vortices,  so  generated  move  out  from  the  wall 
as  they  are  convected  downstream.  These  small  disturbances  ultimately 
become  part  of  the  structure  of  the  turbulence  in  the  wall  shear  layer 
and  ultimately  become  responsible  for  the  Reynolds'  stresses  of 
turbulent  flow. 

The  idea  fundamental  to  this  theory  is  that  large  scale  disturbances 
which  produce  Reynolds'  stresses  some  distances  downstream  were, 
at  an  earlier  time,  small  disturbances  at  the  edge  of  tne  viscous  sublayer 
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some  distances  upstream.  Once  the  small  disturbances  start  to  grow, 
the  e'fect  of  the  polymer  molecules  may  be  neglected.  All  that  the 
polymer  molecules  do  is  to  slightly  alter  the  energy  balance  of  the 
turbulence  fluctuations  close  to  the  wall.  By  alteration,  the  molecules 
allow  viscous  dissipation  to  destroy  disturbances  which  would  have  had 
sufficient  energy  to  grow  had  the  polymer  molecules  not  been  present. 
By  decreasing  the  number  of  small  disturbances  which  grow  per  mit 
area  and  time  and  move  out  from  the  edge  of  the  viscous  sublayer,  the 
addition  of  polymer  molecules  ultimately  changes  the  structure  of  the 
turbulence  in  the  wall  shear  layer.  This  change  results  in  lower 
Reynolds'  stresses  and  thus  lower  turbulent  momentum  transport. 

The  experimental  investigation  of  the  turbulent  structure  of  drag 
reducing  flows  has  been  the  subject  of  study  for  many  researchers. 
Rodriguez,  Patterson,  and  Zankin  (1971)  as  well  as  Wells,  Harkness, 
and  Meyer  (1968)  measured  turbulence  intensities  using  piezoelectric 
probes.  Seyer  and  Metzner  (1969)  measured  mean  velocity  profiles 
and  the  axial  and  radial  turbulence  intensities  for  a  drag  reducing  flow 
using  an  air  bubble  visualization  technique.  Shaver  and  Merrill  (1959) 
used  a  dye  injection  technique  to  qualitatively  measure  the  effect  of 
polymer  addition  on  turbulence.  Virk  et  al  (1967)  and  Patterson  and 
Zankin  (1967)  measured  mean  velocity  profiles  and  axial  intensities 
with  hot  film  anemometers. 
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Nicedemo,  Acierno,  and  Astarita  (1969)  and  Elata,  Lehrer,  and 
Kahanovitz  (1966)  measured  velocity  profiles  using  Pitot  tube  probes. 
More  recently  Rudd  (1972)  measured  the  axial  mean  velocity  profiles 
and  turbulence  intensities  in  a  pipe  with  square  cross-section  using  a 
laser  anemometer  technique.  Using  a  similar  technique,  Chung  and 
Graebel  (1972)  made  these  measurements  in  a  pipe  with  round  cross- 
section,  as  well  as  energy  spectra  at  various  radial  stations. 

The  results  of  many  of  these  experiments  are  inconclusive.  Much 
of  the  data  shows  a  large  degree  of  scatter.  Seyer  and  Metzner  (1969) 
have  observed  significant  differences  in  the  axial  turbulence  intensities 
of  the  polymer  solution  as  compared  to  a  water  flow.  However,  these 
differences  are  of  the  same  order  of  magnitude  as  their  confidence  band. 
Thus,  it  is  not  clear  whether  these  noted  differences  are  real  or  only 
apparent. 

Many  of  these  experimental  results  are  also  questionable  due  to  the 
measuring  techniques  used.  Astarita  and  Nicodemo  (1966)  among  others 
have  questioned  the  accuracy  of  measurements  made  using  hot  film 
anemometers  and  Pitot  tube  probes.  Pitot  tube  errors  are  attributed  to 
the  fact  that:  (a)  the  influence  of  the  normal  stress  is  significant; 

(b)  the  time  average  of  the  fluctuating  stresses  is  not  simply  related 
to  the  time  average  of  the  velocity  fluctuations;  (c)  the  boundary-layer 
hickness  on  the  pressure  probe  may  be  large  compared  to  the  probe 
size;  (d)  a  large  size  Pitot  tube  is  needed  to  minimize  calibration 
errors  (Metzner  and  Astarita,  (1967)  ). 
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Hot -film  anemometers  have  different  heat  transfer  characteristics 
in  pi  lymer  solutions  than  in  water  and  therefore  cannot  be  calibrated 
using  the  usual  heat  transfer  laws.  They  further  suffer  from  the  problem 
of  a  time  dependant  calibration  due  to  the  buildup  of  polymer  molecules 
near  the  stagnation  region  of  the  probe  (Lumley,  1969). 

It  is  not  surprising,  then,  that  these  investigators  report  contra¬ 
dictory  results.  For  example,  Virk  et  al  (1967)  measured  higher  axial 
turbulence  intensities  for  the  polymer  solution  across  a  pipe  flow,  as 
compared  to  solvent  flow,  whereas  Seyer  and  Metzner  (1969)  measured 
lower  intensities  in  the  same  region  of  a  pipe  flow.  On  the  other  hand, 

Wells,  Harkness,  and  Meyer  (1968)  conclude  from  their  results  that 

the  axial  turbulence  intensity  is  unchanged  for  the  polymer  solution  \ 

» 

l 

across  the  pipe  cross-section. 

Most  authors  agree  that  there  is  some  sort  of  interaction  between  1 

the  polymer  molecules  and  the  turbulent  eddies  such  that  the  turbulent  : 

flow  structure  is  altered.  It  is  further  generally  agreed  that  the  pre¬ 
dominant  alteration  occurs  near  the  solid  boundary,  that  is,  in  the 
vicinity  of  the  viscous  sublayer.  There  is,  however,  little  data 
documenting  the  turbulence  processes  close  to  the  wall. 

The  experimental  picture  would  be  bleak  indeed,  if  it  were  not  for 

s 

the  development  of  the  laser-Doppler  anemometer.  In  the  laser- 
Doppler  anemometer  one  projects  a  beam  of  laser  light  through  a  flow 
into  which  minute  particles  have  been  injected.  These  particles  serve 
to  scatter  the  incident  light,  inducing  a  Doppler  frequency  shift  in  the 
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scattered  light  which  is  proportional  to  the  fluid  velocity.  Further 
details  of  this  process  will  be  discussed  later.  Using  this  technique, 

Rudd  (1972)  has  overcome  many  of  the  experimental  problems  inherent 
in  measuring  polymer  flows  and  had  measured  the  axial  mean  velocity 
and  turbulence  intensity  profiles  well  into  the  viscous  sublayer.  Thus 
he  was  able  to  verify  the  growth  of  the  laminar  sublayer  which  was 
predicted  by  Elata,  Lehrer,  and  Kahanovitz  (1966).  (The  growth  of  the 
laminar  sublayer  in  polymer  flows  was  also  observed  indirectly  by 
Virk  (1971a)  ).  This  growth  is  explained  to  be  the  consequence  of 
damping  of  high  wave  number  turbulence  by  the  polymer  molecules. 

The  laser  anemometer  measuring  technique  shows  a  great  deal 
of  promise  as  it  has  many  inherent  advantages  over  traditional  fluid 
flow  measuring  devices  such  as  hot-wire  and  hot  film  anemometers  and 
Pitot  tube  probes.  It  is  more  accurate  over  a  wide  range  of  velocities; 
its  response  is  linear  to  any  component  of  velocity  independent  of  other 
components;  its'  absolute  characteristics,  determined  by  optical 
configuration  and  the  frequency  of  the  laser  light  used,  eliminate  the 
need  for  calibration  under  operating  conditions;  it  has  a  faster  response 
to  turbulent  fluctuations;  and  it  offers  greater  spatial  resolution  which 
permits  a  highly  localized  probe  of  velocity  without  perturbation  of  the 
flow  (Jackson  and  Paul,  1971). 

Since  its  first  development  and  use  by  Yeh  and  Cummins  (1964),  the 
laser-Doppler  anemometer  (hereafter  abbreviated  as  LDA)  has  experienced 
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rapid  development.  Foreman  et  al  (1965)  devised  a  practical  means 
of  detecting  the  time  varying  Doppler  frequency  shift. 

Goldstein  and  Kried  (1967)  used  the  LDA  to  measure  flow  develop¬ 
ment.  Their  optical  configuration  offers  some  decided  advantages  over 
that  of  Foreman  et  al. 

Rolfe  et  al  (1968)  developed  design  criteria  for  the  LDA  and 
constructed  a  LDA  capable  of  measuring  turbulence  properties  of  a  flow 
in  all  three  coordinate  directions.  One  outgrowth  of  Rolfe's  work  has 
been  the  design  and  manufacture  of  a  wideband  frequency  tracker  by 
Fridman,  Kinnard,  and  Meister  (1969),  which,  as  with  Foreman's  device, 
solves  the  knotty  problem  of  extracting  the  real-time  velocity  information 
from  the  time  varying  output  of  the  LDA. 

Because  of  its  many  advantages  and  its  proven  success  in  overcoming 
the  problems  inherent  in  making  measurements  in  polymer  flows,  a 
LDA  was  used  in  the  present  work.  Using  this  technique,  turbulence 
measurements  of  drag  reducing  flows  were  made  in  a  two  inch  diameter 
pipe.  The  goal  of  the  present  work  was  to  extend  our  knowledge  of  the 
turbulent  structure  of  drag  reducing  flows  in  two  ways:  first,  by 
documenting  turbulence  parameters  not  measured  previously,  such  as 
the  Reynolds'  stress,  autocorrelation  function,  radial  energy  spectra, 
and  the  instantaneous  product  of  the  velocity  fluctuations.  (This  last 
measurement  was  made  in  order  to  learn  more  about  the  production  c: 
turbulent  energy.)  Secondly,  measurements  of  all  the  turbulence 
parameters  documented  (mean  velocity,  axial  and  radial  intensities, 
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Reynolds*  stress,  energy  spectra,  the  autocorrelation  function,  and 
the  product  of  the  velocity  fluctuations)  were  made  in  the  region  close 
to  the  wall.  As  mentioned  previously,  most  authors  feel  that  the  poly¬ 
mer-turbulence  interaction  that  results  in  drag  reduction  occurs  in  this 
region.  It  is  hoped  that  the  results  of  these  measurements  will  add 
new  insight  into  the  basic  nature  of  the  drag  reduction  phenomenon. 


Chapter  It 


THE  PRINCIPLE  OF  OPERATION 
OF  THE  LASER-DOPPLER  ANEMOMETER 

The  principle  of  operation  of  a  laser-Doppler  anemometer  can  be 
briefly  described  as  follows:  an  incident  beam  of  laser  light  (in  the 
direction  designated  by  the  unit  vector  kQ  in  Figure  2. 1),  is  focused 
at  point  0  in  a  fluid  stream  whose  velocity  at  0  is  V.  The  incident  beam 
is  scattered  due  to  the  presence  of  minute  particles  injected  into  the 
flow.  The  particles  move  with  nearly  the  velocity  V  (Huffaker,  1970). 

A  portion  of  the  scattered  front,  a  beam  denoted  by  k  in  Figure  2.  1, 
has  a  Doppler- shifted  frequency  f  =  f  +  f^,  where  fg  is  the  frequency 
of  the  scattered  light,  f  is  the  frequency  of  the  incident  light,  and  f^ 
is  the  Doppler  shift  due  to  scattering.  The  Doppler-shifted  frequency, 
fd,  is  directly  proportional  to  the  velocity  of  the  scattering  particles 
at  point  0  (the  exact  relationship  is  developed  below),  thus  fd  is  a  measure 
of  the  stream  velocity  Vat  point  0.  In  order  to  detect  fd,  a  reference 
beam  of  the  initial  frequency  is  aligned  so  as  to  coincide  with  the  scattered 
beam  and  is  focused  at  point  0.  The  scattered  and  reference  beams  are 
optically  mixed  at  point  0.  The  resulting  beam  is  optically  heterodyned 
on  the  photo  cathode  surface  of  a  photomultiplier  tube  to  produce  a  beat 
signal  with  a  frequency  of  the  Doppler  shift,  f  .  Thus  the  photomultiplier 
(P.  M.)  output  signal  is  a  voltage  whose  frequency  is  the  Doppler-shifted 
frequency. 
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Figure  2. 1 

Optical  Geometry  for  a  Laser-Doppler  Anemometer 


The  relationship  between  the  Doppler  frequency  shift,  f^,  and  the 
velocity  Vcan  be  found  in  the  following  manner:  the  incident  frequency, 
f',  as  seen  by  an  observer  moving  with  the  flow,  is 

V  •  k 

(2.  1)  f*  =  fQ  (1  -  , 

where  C  is  the  velocity  of  light  in  the  fluid,  and  k  Q  and  k  are  the  unit 

vectors  in  the  direction  of  the  incident  and  scattered  beams  respectively. 

The  frequency  f  of  light  scattered  by  particles  moving  with  the 
s 

flow  and  received  by  a  stationary  detector,  say  the  P.M. ,  is 

V  •  k 

(2.2)  f.  =  f’(l+  — g— 1 ■)  . 


It  follows  that 

V  •  k  V  •  k 

<2*3)  fs  =  V1-  ZncZ2)(1  +  • 


The  Doppler  frequency  shift,  f^,  then  becomes 


f.  =  f  -  f 
d  so 


(2.4)  f .  =  f  [—  •  (k  -  k  )  -  ( 
d  o  L  r  — s  —  o' 


V  •  k  V  •  k 

)]• 


.  v  . 

which  to  order  |  |  is 


(2.  5) 


fd  -  -t  ?  •  fes-io>- 


The  optical  arrangement  used  by  Foreman  et  al  (1965)  is  shown 


in  Figure  2.  2  and  that  of  Goldstein  and  Kreid  (1967)  in  Figure  2.3. 
Because  of  its  relation  ease  of  alignment,  a  two-component  variation 
of  Goldstein's  arrangement  was  used  in  this  work  (see  Figure  3. 1). 
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Mirror 


Figure  2.  2 

Optical  Arrangement  of  Foreman  et  al  (1S66) 
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Flow 


Figure  2.3 

Optical  Arrangement  of  Goldstein  and  Kreid  (1967) 
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j  Flow  Direction 


Figure  2.4 

Optical  Geometry  for  the  Two-Component  Laser-Doppler  Anemometer 
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The  relationship  between  the  components  of  V  and  the  Doppler- 
shifted  frequencies  for  our  two-component  optical  geometry  can  be 
determined  by  considering  equation  (2.  5)  and  Figure  2.4.  The  flow 
velocity  V  in  terms  of  the  cylindrical  coordinates  defined  in  Figure 
2.4  is 


(2.6) 


I  -  Vrer  +  Vzez  +  V,efl. 


The  direction  of  the  incident  beam  is  defined  by  the  unit  vector  k 

— o 


(2.7) 


k  -  -  e 
-o  — r 


The  direction  of  a  ray  of  laser  light  scattered  through  an  angle  9  above 
the  incident  beam  is  defined  by  the  unit  vector  k  ^  ,  where 


(2.8) 


-si  =  '  cos0£r  "  sine  £z  • 


Likewise  the  direction  of  a  ray  of  laser  light  scattered  through  an 
angle  9  below  the  incident  beam  is  defined  by  the  unit  vector  k  g2  , 


where 


(2.9) 


k  „  =  -  cosge  +  singe 
—  s2  -r  — z 


The  Doppler-shifted  frequency,  fdp  in  the  kg^  direction  is,  using 
equation  (2.  5) 


f  V 

o- 


(k  _  „  -  k  ) 

—  si  —  o'  ’ 


_  f0 


C“<vr£r+  vz-z  +  ve-e)  '  f(1"cos0)£r-sin0£z] » 


Chapter  in 

THE  OPTICAL  SYSTEM 

The  two  most  important  parameters  governing  the  performance 
of  an  experimental  system  are  the  signal  to  noise  ratio,  SNR  ,  and 
the  frequency  response.  In  order  to  obtain  a  high  over-all  SNR  in  the 
laser-Doppler  anemometer,  the  optical  components  must  be  carefully 
aligned  in  accordance  with  some  stringent  optical  requirements. 

(These  requirements  will  be  developed  in  the  next  chapter.)  The  key 
to  the  operation  of  the  LDA  as  a  useful  experimental  tool  is  an  optical 
bench  that  will  permit  easy  alignment  and  adjustment  of  its  components. 
The  optical  bench  used  in  this  work  was  designed  and  built  at  The 
University  of  Michigan  in  order  to  satisfy  the  optical  requirements  and 
also  to  offer  ease  of  alignment  and  adjustment. 

1.  Optical  Arrangement 

The  optical  arrangement  used  in  this  experiment  (Figure  3.  1)  was 
a  variant  of  the  one  used  by  Goldstein  and  Kreid  (1967).  It  can  be  briefly 
described  as  follows:  the  output  beam  of  a  Helium-Neon  laser  (beam 
Bq  in  Figure  3. 1)  is  split  at  beam  splitter  BS^  (10%  reflection)  into 
beams  Bq  and  B  Beam  Bp  a  reference  beam,  is  reflected  at  mirror 
M,  and  focused  at  point  0  by  means  of  lens  L.,  where  it  forms  an  angle 
9  with  respect  to  beam  BQ.  After  passing  through  point  0,  beam  B ^ 
passes  through  the  front  aperture  Ajp  the  front  receiving  lens  L^, 
pinhole  P^,  the  rear  receiving  lens  Lf  y  and  the  rear  aperture  Arl, 
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before  striking  photomultiplier  P.M.  y  The  reference  beam  is 
optirally  heterodyned  at  the  P.M.  photocathode  surface  with  the  portion 
of  the  scattered  beam  front  (from  beam  Bq)  allowed  through  apertures 

and  and  pinhole  P  ^  to  produce  a  signal  at  the  P.  M.  whose 
frequency  is  the  Doppler -shifted  frequency.  The  output  signal  produced 
by  the  optics  and  P.  M.  along  beam  B  <  is  termed  the  output  of  channel 
1. 

Similarly  beam  B2,  a  reference  beam,  is  formed  by  beam  splitter 
BS2,  (25%  reflection  ),  reflected  at  mirror  M2  and  focused  by  lens 
L2  at  point  0  where  it  forms  an  angle  e  with  beam  B^.  Beam  B^  then 
passes  through  aperture  Aj2  lens  Lp,  pinhole  P2,  lens  Lr2,  and  strikes 
P.  M.  2  after  passing  through  aperture  Ar2>  Beam  B2  is  heterodyned 
at  the  P.M.  with  a  portion  of  the  scattered  wave  to  produce  another 
Doppler-shifted  signal.  The  output  signal  produced  by  the  optics  and 
P.M.  along  B2  is  termed  the  output  of  channel  2.  After  the  beam 
splitter  BS2,  beam  BQ,  the  scattering  beam,  passes  through  lens  LQ 
which  focuses  it  at  point  0  where  it  is  scattered.  The  trisection  of 
beams  BQ,  B^,  and  B2  at  point  0  forms  the  measuring  volume. 

2.  Optical  Geometry 

The  optical  geometry  can  be  described  as  follows:  Beam  Bq  was 
aligned  so  as  to  be  perpendicular  to  the  outside  faces  and  F2  of  the 
test  section  and  to  pass  through  the  axis  of  flow,  D.  Beams  Bj  and  B2 
were  placed  at  equal  angles  B  on  either  side  of  beam  B^  and  lay  in  the 
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plane  formed  by  the  axis  D  and  beam  Bq  and  intersected  at  point  0.  This 
choice  of  optical  geometry  provided  an  easy  method  of  obtaining  the 
axial  and  radial  components  of  velocity.  (In  the  same  way  as  one  obtains 
these  quantities  from  the  output  of  a  X-wire  in  hot-wire  anemometry. ) 
This  configuration  further  minimized  optical  distortion  due  to  the 
presence  of  the  test  section.  This  geometry,  a  two-component  variation 
of  the  one  used  by  Goldstein  and  Kreid  (1967),  offered  an  additional 
advantage  over  the  geometries  used  by  Foreman  et  al  (1966)  and  Rolfe 
et  al  (1968)  in  that  the  optical  mixing  of  the  scattered  and  reference 
beams  in  the  present  system  took  place  at  the  measuring  volume  rather 
than  at  the  P.M.  Therefore,  beam  alignment  became  less  critical. 

3.  Optical  Comp  ments 

All  of  the  optical  components,  lenses,  mirrors,  photo  multiplier 
tube  holders,  and  the  laser,  were  mounted  on  a  vertical  rigid  reference 
plate,  14  inches  high  by  66  inches  long.  This  plate  was  formed  from 
1/8  inch  thick  steel  plate  to  the  above  dimensions  with  1  3/4  inch 
flanges  top  and  bottom.  The  rigidity  of  the  plate  was  increased  by 
welding  a  2  l/2  x  2  l/2  x  66  inch  angle  iron  along  its  length,  7  1/2 
inches  up  from  its  base.  The  plate  was  mounted  to  a  300  pound  21  by  31 
inch  granite  surface  plate  by  means  of  four  triangular  braces  and  eight 
screws.  The  resulting  reference  plate-surface  plate  system  was  thus 
very  rigid. 
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In  order  to  obtain  measurements  across  the  test  section,  point  0, 
the  measured  point,  (Figure  3.  1)  must  be  moved  relative  to  the  test 
section.  This  was  accomplished  by  a  set  of  six  rollers,  placed  between 
the  granite  table  and  the  bench.  Translation  of  the  measuring  point 
was  accomplished  by  means  of  a  spring  loaded  screw,  one  end  mounted 
to  the  bench,  the  other  to  the  granite  table.  (See  Figure  3.  1)  This 
allowed  the  entire  optical  system  (and  therefore  point  0)  to  be  moved 
as  a  unit.  The  transverse  position  of  the  measuring  point  0  in  the  test 
section  was  determined  by  first  locating  it  visually  (using  a  telescope) 
at  the  pipe  center.  Subsequent  positions  were  determined  by  reading 
a  dial  indicator  placed  against  the  granite  table  in  the  line  of  motion. 

Due  to  the  refraction  of  the  beams  by  the  test  section,  the  distance 
point  0  moved  in  the  test  section  was  4/3  of  the  dial  indicator  movement. 
(Details  of  the  refractive  effect  of  the  test  section  on  beam  geometry 
can  be  found  in  Appendix  A. ) 

The  bench  mentioned  above  consisted  of  1  1/2  inch  plywood  mounted 
on  Unistrut  supports  which  were  in  turn  mounted  on  a  cement  block 
base.  Asphalt  saturated  felt  was  placed  between  the  cement  blocks  to 
isolate  the  optical  bench  from  building  vibrations. 

The  test  section  was  machined  from  a  solid  piece  of  3.  125  by  3.  125 
by  4  inch  Plexiglas  (Index  oi  refraction  1.  5).  A  two  inch  hole  was 
machined  through  the  center  of  the  square  cross  section.  It  was  recessed 
at  either  end  to  fit  smoothly  into  the  two  inch  internal  diameter  brass 
pipe  to  which  it  was  attached  by  flanges.  The  optical  surfaces  of  the 
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test  section  were  polished  to  a  scratch-free  finish  by  means  of 
successive  applications  of  Bon-Ami  and  Cadco  Anti-Static  polish. 

The  inner  surface  of  the  test  section  was  cleaned  periodically  of 
a  build-up  of  dirt,  latex,  and  polymer  to  avoid  impairment  of  its  optical 
qualities.  This  cleaning  was  accomplished  by  removing  the  lower 
flange  from  the  test  section  so  that  the  inner  surface  could  be  cleaned 
with  Cadco  Anti-Static  polish.  The  polish  was  applied  with  lint-free 
cotton  balls  on  the  end  of  a  dowel  rod. 

The  laser  (Spectra-Physics  Model  120,  Helium-Neon  laser,  wave¬ 
length  6328  angstroms)  was  mounted  on  an  I-beam  which  was  attached 
rigidly  to  the  vertical  steel  reference  plate.  The  mounting  holes  in  the 
I-beam  were  slotted  to  allow  adjustment  of  the  laser  to  pass  through 
the  center  of  the  test  section  and  also  to  be  perpendicular  to  its  outer 
surfaces. 

The  front  surface  mirrors  used  in  this  experiment  were  mounted 
in  he  mi-spherical  holders.  These  holders  were  attached  firmly  to  the 
mirror  mounts  with  magnets.  The  mirror  mounts  were  then  attached 
to  the  reference  plate  by  slotted  bars  and  held  in  position  by  knurled 
tightening  nuts.  (See  Figure  3. 2)  The  slotted  bars  allowed  for  gross 
translation  in  two  dimensions.  Tne  mirror  mount  could  be  fine  adjusted 
in  one  direction  by  means  of  a  screw  operated,  spring  loaded  transla¬ 
tion  stage.  This,  along  with  the  two  degrees  of  rotational  freedom 
allowed  by  the  hemi  spherical  mirror  holders,  permitted  easy  mirror 
adjustment. 


The  beam  splitter  holders  were  mounted  perpendicular  to  the 
ref<  rence  plate  by  a  sleeve  attached  to  the  plate.  (See  Figure  3.2.) 

The  angle  between  beam  splitter  holder  and  beam  could  be  adjusted  by 
rotating  the  holder  pin  in  the  sleeve.  It  was  held  in  place  by  a  knurled 
tightening  nut.  Fine  angular  adjustment  of  the  beam  splitter  was 
accomplished  by  adjusting  the  set  screws  located  in  the  holder. 

The  lens  holders  were  mounted  on  °  slide  array  which  allowed 
for  gross  translation  in  the  three  coordinate  directions.  The  slide  array, 
which  was  mounted  to  the  reference  plate,  consisted  of  a  series  of 
spring-loaded  blocks  mounted  on  drill  rod  as  shown  in  Figure  3. 2. 

The  lens  holders  were  mounted  on  several  blocks  of  the  same  design. 

The  lens  holders  themselves  could  be  fine  adjusted  in  two  coordinate 
directions  by  screw  operated,  spring-loaded  translation  stages  (similar 
to  those  used  in  the  mirror  mounts).  A  yoke  mechanism  provided  an 
additional  two  degrees  of  rotational  freedom.  These  adjustments 
were  found  sufficient  for  accurate  lens  placement  and  alignment. 

The  photo multiplier  tube  holder  mounts  were  attached  to  the  refer¬ 
ence  plate  on  one  side  and  to  the  supporting  table,  with  vertical  rods, 
on  the  other  side.  (See  Figure  3.  3)  This  mechanism  allowed  for  gross 
translation  of  the  P.M.  holders  in  the  three  coordinate  directions 
utilizing  the  same  type  of  slide  array  design  used  in  the  lens  system. 
Yokes  mounted  on  the  front  of  the  slide  array  allowed  for  rotation  in 
tvo  directions  about  the  leading  aperture,  A^j,  in  the  P.M.  holder. 

A  yoke  mechanism  mounted  on  the  rear  slide  array  likewise  provided 
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the  rotation  about  two  coordinate  axis  needed  for  gross  alignment 
of  the  P.M.  holders.  Finally,  spring-loaded  screw  adjustable  adjust¬ 
ment  frames  mounted  on  the  rear  slides  provided  fine  angular  adjustment 
of  the  P.M.  holder  about  the  leading  aperture. 

The  P.M.  holder  (see  Figure  3.3)  housed  the  P.M.  tube  as  well 
as  the  apertures  and  Ar  and  the  lens-pinhole  combination  of  Lp 
and  P.  The  apertures,  which  were  adjustable,  prevented  extraneous 
light  from  reaching  the  P.M.  tube.  They  also  provided  an  easy 
reference  for  aligning  the  P.M.  holder  along  the  beam  axis  as  well  as 
providing  a  means  to  block -off  the  unwanted  portion  of  the  reference 
beam.  The  pinhole-lens  combination  (suggested  by  F.  Danel)  was 
designed  so  that  the  image  of  the  measuring  point  0  was  projected  on 
the  P.M.  tube.  The  pinhole  diameter  was  the  same  size  as  the  measur¬ 
ing  volume  cross-section.  Therefore,  any  light  extraneous  to  the 
measuring  volume  was  cut-off  by  the  pinhole.  Thus  the  P.M.  "saw" 
only  the  measuring  volume  and  the  error  due  to  sensing  velocities 
from  other  points  was  minimized.  This  became  especially  important 
when  making  measurements  close  to  reflecting  boundaries.  This  design 
was  instituted  when  it  was  fc  ind  that  the  heterodyned  signal  was  too 
"noisy"  when  the  measuring  volume  was  close  to  the  pipe  wall.  The 
pinhole -lens  combination  was  mounted  as  a  unit  and  could  be  moved 
along  the  P.  M.  holder  axis  by  means  of  a  spring-loaded  adjusting 
screw.  This,  along  with  the  adjustable  bushing  adapter  (see  below) 
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allowed  the  pinhole-lens  combination  to  be  adjusted  to  produce  the 
be  s’,  signal. 

The  main  body  of  the  P.M.  holder  was  made  from  two  sizes  of 
copper  tubing  which  also  served  as  an  effective  electrostatic  and 
magnetic  shield  (see  Figure  3.3). 

The  P.  M.  holder  was  completely  light-sealed  and  had  a  cover 
which  fit  over  the  front  aperture  to  prevent  background  light  from 
entering  and  fatiguing  the  P.M.  tube  when  the  LDA  was  not  in 
operation.  The  inside  surfaces  of  the  holder  were  painted  ultra-flat 
black  to  absorb  any  stray  light. 

The  P.M.  tube  socket  was  contained  in  a  housing  along  with  the 

voltage-dividing  circuitry  for  the  P.M.  tube.  The  housing  had  B.N.C. 

connections  for  input  from  the  power  supply  and  output  to  the  readout 

unit.  (This  output  was  the  heterodyned  signal. )  The  housing  slid  into 

a  bushing  adapter  which  in  turn  slid  into  the  P.  M.  holder.  The  bushing 

adapter  could  be  moved  along  the  P.M.  holder  axis  so  that  the  P.M. 

front-surface-to-lens  L  distance  could  be  varied.  The  P.M.  tube 

r 

and  circuit  assembly  were  held  in  place  in  the  P.M.  holder  with  set 
screws.  This  allowed  for  easy  removal  of  the  assembly  to  allow 
adjustment  of  the  P.  M.  holder  during  alignment. 

All  reflecting  surfaces  of  the  optical  bench  were  painted  ultra-flat 
black  to  minimize  the  incidence  of  stray  light.  Furthermore,  the  entire 
optical  bench  was  housed  in  a  light-sealed  room  to  isolate  it  from  the 
main  laboratory  area. 


Chapter  IV 


SIGNAL-TO-NOISE  CONSIDERATIONS 
FOR  THE  LASER-DOPPLER  ANEMOMETER 

The  laser-Doppler  anemometer  uses  the  proportionality  of  the 
Doppler -shifted  frequency  of  the  scattered  wave  to  the  fluid  velocity 
at  the  measuring  point  to  make  flow  measurements.  The  anemometer 
used  in  this  experiment,  described  previously,  mixed  the  Doppler- 
shifted  scattered  wave  with  a  reference  wave  of  the  initial  unshifttd 
frequency  at  the  measuring  point,  0.  The  resulting  wave  was  sensed 
by  a  photomultiplier  tube  (RCA  7326).  The  photomultiplier  can  be 
operated  in  two  modes;  as  a  detector  (low  level)  of  electromagnetic 
waves  or  as  a  mixer  of  electromagnetic  waves,  as  in  the  present 
experiment.  As  used  in  the  current  work,  the  photomultiplier  acted 
as  a  superheterodyne  receiver,  hi  this  mode  of  operation  of  the  P.M. , 
the  frequency  of  the  incoming  Doppler-shifted  signal  was  mixed  with 
that  of  the  unshifted  reference  signal  (local  oscillator)  to  a  new,  lower, 
frequency.  This  signal  was  then  amplified  before  further  processing 
was  done  by  instrumentation  subsequent  to  the  P.M.  The  P.M. 
possesses  advantages  over  other  sensing  devices  in  that  it  offers 
an  essentially  noiseless  post  detection  gain.  Furthermore,  the  band¬ 
width  of  the  P.M.  allows  measurement  of  velocities  into  the  supersonic 
range,  far  beyond  the  needs  of  this  experiment. 
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1.  Noise  in  the  Pt  M.  Tube 

Noise  in  the  P,M.  tube  is  generated  by  random  electron  motion. 

It's  spectrum  is  flat,  i.e. ,  the  noise  is  white.  The  three  dominant 
noise  sources  are:  (1)  Dark  current.  When  the  P.  M.  is  operated  in 
the  absence  of  incident  light,  an  anode  current  flows  which  sets  a  limit 
to  the  lowest  intensity  of  light  that  can  be  measured.  It  is  caused  mainly 
by  thermionic  current  emitted  from  the  photosensitive  coating.  (2)  When 
light  impinges  on  the  P.M.  tube,  the  output  noise  increases  significantly 
because  of  statistical  fluctuations  of  the  photoemission  current  modified 
by  the  gain  variation  of  the  multiplier.  This  is  called  shot  noise  or 
"noi..e-in-signal"  since  it  is  a  function  of  the  input  signal.  (3)  Shot 
noise  contributed  by  extraneous  light  (light  from  sources  from  other 
than  the  measuring  volume)  (Rolfe  et  al,  1968). 

The  noise  from  background  light  can  be  minimized  by  placing 
apertures  in  front  of  the  P.M.  tube  to  block  off  this  light.  The  dark 
current  is  negligible  in  comparison  to  the  current  due  to  the  incident 
laser  beam  light.  Therefore,  the  dominant  noise  source  which  limits 
the  sensitivity  of  the  LDA  is  the  noise  in  the  signal. 

2.  Noise  in  the  Signal 

The  two  most  important  parameters  governing  the  sensitivity  of 
an  instrument  are  its  signal-to-noise  ratio  and  its  bandwidth  (frequency 
response).  In  this  section  an  expression  for  the  signal-to-noise  ratio 
(SNR)  for  the  P.M.  tube  is  developed. 
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The  first  step  in  this  development  is  to  calculate  the  P.M.  anode 
current. 

Consider  two  light  beams  incident  on  the  photocathode  surface: 

E  sinw  t  due  to  the  scattered  radiation  and  E  sin  w  t  due  to  the 
s  s  r  r 

reference  radiation,  where  E  denotes  the  electric  field  amplitude  and 
w  is  the  radian  frequency  of  the  light.  These  waves  add  to  become  the 
total  incident  radiation  at  the  cathode  whose  electric  field  is 

(4.1)  E,  =  E  sinw  t  +  E  sinw  t  , 

v  '  t  s  s  r  r  * 

assuming  the  beams  have  the  same  linear  polarization  direction. 

Since  the  P.M.  is  a  square  law  detector,  the  anode  current  out¬ 
put  is  proportional  to  the  power  of  the  incident  light  beam 


(4.2)  I.  =  S  HE  2 sin2 w  t  +  E  2sin2w  t  +  2E  E  sinw  tsinw  t] 

'  ta  Ls  s  r  r  rsr  sJ 


where  the  current  amplification  S  is  assumed  constant.  This  can  be 


rewritten  as 

.  t  cos2w  t  cos2w  t 

(4.3)  I  =  S  {E2(5  -  —j - )  *  E  2  (i - j-i-) 


+  EjEg  [cos  (ws  -  wr)  t  -  cos  (wg  +  wr)  t  ]  }  . 


The  terms  containing  the  frequencies  2wg,  2wf,  and  (wg  +  are  not 
important  since  they  fall  outside  the  bandwidth  of  the  P.M.  Therefore, 
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E  +  E 

(4.4)  Ita  =  S  { — S  2  r  +  EsErcos(ws  -  wr)  t} 
or 

(4.5)  1^  =  Sk  {Pg  +  Pr  +  cosAwt  }  , 


where  P  is  the  light  power  impinging  on  the  photocathode,  Aw  is 

w  -  w  ,  and  S.  is  the  cathode  radiant  sensitivity.  This  is  the  total 
s  r  K 

current  emitted  by  the  cathode  in  response  to  the  laser  light. 

The  noise  current,  i^,  is  given  by  the  shot  noise  equation 


(4.6)  i^2  =  2eAfIta  , 

(4.7)  i^  =  2eAfSk[Ps+Pr  +  2^PTP^cosAwt]  , 

where  e  is  the  electron  charge.  This  represents  the  noise  current 
spread  randomly  over  the  frequency  range  of  bandwidth  Af.  This 
noise  can  be  measured  by  measuring  its  power, 

<ik2>R  , 

2 

where  ik  is  the  mean  square  noise  current  and  R  is  the  resistive 
component  of  the  circuit  impedenee.  Then 


(4.8) 


..  2 


<*k  >  “  2eAfVps  +  IV  ’ 


2 

where  we  have  averaged  i.  to  get  the  mean  square  noise  current. 

K 

The  a.c.  current  at  the  cathode  is,  as  we  have  seen 


I  =  2S,  -5P~P"  cos  Awt  . 
ac  k v  s  r 


(4.9) 
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The  mean  square  of  the  signal  current  is 
(4.10) 

and  the  cathode  SNR  is 


<1  2  >  =  2S 2 P  P  , 
ac  k  s  r  ’ 


(4.11) 


<ij> 


SkP_P 


s  r 


e  Af  (Pg  +  Pr) 


The  signal-to-noise  ratio  at  the  anode  becomes 


(4.12) 


where 


/  S  a  _  1  /  S  v  __ 

'  X.  'a  '  XT  '/»  — 


S^P  P 
k  s  r 


Ki'a  "  ¥  ^  N'c  “  e  Afk  (P  +Pl  » 

v  s  r' 


a 

a  -1  » 


a  being  the  average  gain  per  stage  of  the  multiplier  chain. 
Rev.  riting  equation  (4.12)  we  have, 

(4.13) 


(S-)  -  S  p 

^N'a  e  Af  k  s 


1  + 


We  see  that  making  P  >>  P  will  maximize  th*  SNR  until,  in  the 

r  s 

limit  as  P  /'P  approaches  infinity, 

JT  S 


(4.14) 


(S_) 

'  N;a 


.  p 

e  Afk  s 


For  this  limiting  case  the  SNR  is  seen  to  be  independent  of  Pr . 


However,  P  cannot  be  increased  without  limit.  P  is  restricted 
*  r  r 

in  value  such  that  1^  (see  equation  (4. 5))  does  not  exceed  its  maximum 
value  as  defined  by  manufacturer  specifications. 
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3,  Optical  Loss  Coefficients 

The  preceeding  derivation  of  the  SNR  did  not  take  into  account 
losses  due  to  the  optical  system.  These  can  be  included  by  modifying 
equation  (4.13)  by  introducing  loss  coefficients  in  the  following  manner, 


(4.15) 


CdCtCc  t  e*ATk 


P 

s 


) 


where 

Cd  =  P.M.  alignment  coefficient, 

Ct  =  Transmission  path  coherence  loss  factor, 
and 

=  Source  coherence  length  factor. 

Since  each  of  these  effects  tends  to  degrade  the  SNR,  their  coefficients 
must  be  less  than  or  equal  to  one. 

The  source  coherence  length,  L,  a  parameter  that  defines  the 
coherence  properties  of  a  light  source,  is  defined  as  that  distance 
between  two  wave  packets  emanating  from  the  same  source  such  that 
their  phase  difference  is  A 0„ 

Born  and  Wolf  (1964)  present  an  expression  for  source 
coherence  length  for  A0  =  2i r, 


(4. 16) 


L  =  C/  A  v  , 


where  A  v  is  the  frequency  spread  and  C  the  free  space  velocity  of 
light. 
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The  source  coherence  length  for  lasers  is  on  the  order  of  (10) 

feet.  Therefore,  for  the  present  experiment,  the  source  coherence 

factor,  C  ,  is  essentially  unity, 
c 

In  their  passage  from  laser  to  P.M. ,  the  laser  beams  may  be 
degraded  in  such  a  way  as  to  affect  their  ability  to  produce  coherent 
heterodyne  signals  at  the  P.M.  The  two  causes  of  this  degradation 
are:  (1)  fluctuations  in  effective  path  length  due  to  fluctuations  in 
path  density  and  therefore  refractive  index  and  wavelength;  and  by 
fluctuations  in  geometrical  path  length  with  bent  rays;  (2)  scattering 
from  atmospheric  or  other  sources  outside  the  measuring  volume 
(Rolfe  et  al,  1968).  One  measure  of  these  effects  is  the  transmission 
path  coherence  length.  If  the  difference  between  the  effective  beam 
lengths  of  the  various  beams  used  in  an  experiment  is  less  than  the 
transmission  path  coherence  length,  coherent  heterodyning  will 
be  possible.  Since  the  beam  lengths  of  the  current  work  differed  by 
less  than  4  inches,  the  transmission  path  coherence  length  for  this 
optical  system,  no  coherence  problem  was  encountered. 

Due  to  errors  in  the  optical  mixer  (P.M.)  geometry  and  alignment 
with  respect  to  the  incoming  beams,  the  SNR  can  be  attenuated  signifi¬ 
cantly  by  an  amount  C^. 

In  this  experiment,  alignment  problems  arose  due  to  angular  mis¬ 
match  between  the  incoming  beams  and  the  P.M.  cathode  face  and  also 
due  to  unwanted  beam  overlap  on  the  cathode  face  between  the  scattered 


and  reference  beams. 
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Rolfe  et  al  (1968)  have  shown  that  the  angular  alignment  between 
beams  and  the  P.M.  surface  is  not  critical.  Furthermore,  the  P.M. 
tube  holders  and  their  mounts  used  in  this  work  were  designed  to 
minimize  this  problem. 

That  portion  of  the  reference  beam  that  does  not  intersect  the 
scattered  beam  on  the  P.M.  cathode  face  induces  shot  noise  without 
increasing  the  signal  strength.  Consequently,  this  reference  beam 
overlap  reduces  the  SNR  significantly.  To  avoid  this  problem, 
apertures  and  a  pinhole  were  placed  in  front  of  the  P.M.  to  block-off 
the  unwanted  portion  of  the  reference  beam. 

The  optical  geometry  used  in  the  current  work  offered  some 
distinct  advantages  over  those  used  by  other  investigators.  In  the 
geometries  used  by  Rolfe  et  al  (1969)  and  Foreman  et  al  (1965),  mixing 
of  scattered  and  reference  beams  took  place  on  the  P.M.  tube  photo¬ 
cathode  surface  necessitating  careful  alignment  to  avoid  angular 
mismatch  between  the  beams.  Any  beam  mismatch  would  degrade 
the  SNR.  This  problem  did  not  arise  in  the  current  work  since  the 
scattered  and  reference  beams  mixed  at  their  point  of  intersection 
(point  0  of  Figure  3.1).  Alignment  of  the  P.M.  was  then  effected  by 
merely  aligning  the  P.M.  holder  so  that  the  reference  beam  passed 
through  the  apertures  at  their  smallest  diameter. 

A  further  advantage  of  the  geometry  used  here  is  that  beam  size. 


on  the  P.M. ,  was  immaterial. 
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It  might  be  well  to  mention  here  that  the  laser  anemometer  proposed 
by  Rudd  (1969)  possesses  further  advantages.  This  device*  which  could 
be  termed  an  interference -fringe  anemometer  as  it  does  not  use  the 
Doppler  principle*  has  the  same  mathematical  relationship  between 
detected  frequency  and  flow  velocity  as  the  LDA.  The  velocity  at 
the  measuring  volume  is  detected  as  the  ratio  of  the  distance  between 
interference  fringes  and  the  time  it  takes  for  a  particle  to  pass  between 
them.  The  interference  fringes  occur  at  the  measuring  volume  and 
are  produced  as  follows  (see  Figure  4.1):  Light  from  a  high  power 
laser  is  diverged  by  lens  Lj  to  cover  the  mask  M.  This  light  is 
focused  by  a  single  lens  Lg  to  point  0  in  the  fluid.  Since  both  beams 
pass  through  the  same  lens,  Lg,  they  are  brought  to  focus,  and  there¬ 
fore  overlap,  at  the  same  point  0.  This  system  is  truly  self -aligning. 

An  interference  fringe  pattern  is  formed  at  the  point  0  due  to  the 
interaction  of  the  two  beams  of  equal  intensity  that  intersect  there. 

After  passing  through  the  fluid,  the  beams  are  focused  by  lens  L„  on 
a  photo  detector.  The  output  of  this  detector  is  a  voltage  whose  frequency 
is  related  to  the  velocity  at  the  measuring  volume  in  the  same  way  that 
the  Doppler  frequency  is  in  the  LDA. 

Not  only  does  this  system  provide  ease  of  alignment,  but  the 
frequency  shift  is  independent  of  detector  position.  This  implies  that 
there  is  no  aperture  broadening.  (Aperture  broadening  will  be  dis¬ 
cussed  later  in  this  chapter. ) 


Optical  Arrangement  of  Rudd  (1969) 
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With  all  of  these  advantages,  the  interference -fringe  anemometer 
seems  to  be  a  promising  form  of  laser  anemometer.  It  does  have  the 
disadvantage,  however,  of  requiring  greater  laser  power  with  the 
attendent  safety  problems. 

4.  Instrumental  Broadening 

The  frequency  spectrum  of  the  laser  light  of  a  LDA  will  be 
broadened  (departing  from  purely  monochromatic  light)  due  to  several 
physical  processes: 

(1)  turbulence  phenomena, 

(2)  finite  aperture  size, 

(3)  finite  measuring  volume  size.. 

The  frequency  of  scattered  light  from  a  flowing  fluid  is  sensitive 
to  flow  fluctuations  as  seen  from  the  relationship  developed  in 
Chapter  2, 

(4.17) 

Since  turbulence  is  a  random  phenomenon,  that  is,  it  is  characterized 
by  random  fluctuations  of  velocity  about  some  mean  velocity,  it  is  seen 
from  equation  (4.17)  that  the  frequency  of  the  Doppler -shifted  light 
varies  randomly  about  a  mean  (is  broadened) . 

Anything  that  would  obscure  these  turbulent  fluctuations  (broadening 
from  other  sources)  would  reduce  the  sensitivity  of  the  LDA.  Such 
undesirable  broadening  is  induced  because  of  the  finite  size  of 
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receiving  aperture  and  pinhole  diameters  and  because  of  the  finite 
size  of  the  measuring  volume.  These  effects  will  be  discussed 
below. 

Consider  the  optical  geometry  of  Figure  4.  2.  From  equation 

(4.17) ,  the  Doppler- shifted  frequency,  fd,  is 

2 v 

(4.18)  f,  =  sin  e  . 

d  A 

Therefore 

3*d  2  V 

<4-19>  -irr  =  -F  cos  e  . 

Now,  if  R  is  the  distance  between  the  aperture  and  measuring  point 
and  r  the  aperture  radius,  the  incremental  angular  uncertainty,  A  6  , 
can  be  written  as 

(4.20)  A0  =  . 

The  aperture  broadening  is 

3  fri 

Af  =  A  0  =  A  0  COS0  , 

a  3  0  A 

(4.21)  Afa  =  cos  9  . 

For  this  experiment  the  aperture  broadening  was 

Af  (10) 4  HZ. 
a 

This  was  a  tolerable  level  for  the  present  experiment  since  if)  for  the 

5 

centerline  velocity  was  about  (10)  HZ  at  a  Reynolds’  number  of  50,000. 
However,  this  broadening  must  be  compared  with  the  natural  broadening 
due  to  the  turbulent  fluctuations. 
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The  stream  velocity  V  at  a  point  in  the  flow  can  be  written  as 

(4.22)  V  =  V  +  V' 

where  V  is  the  mean  velocity,  a  constant  in  time,  and  V*  is  the 
variation  of  velocity  about  that  mean. 

Equation  (4. 22)  can  be  written  in  scalar  form  as 


(4. 23)  Vz 

(4. 24)  Vr 

where  z  and  r  denote  the  axial 
and  we  are  not  considering  the 
Recall  that 


=  V  +V’  , 

z  z  ’ 

=  V  +  V'  , 

r  r  ’ 

and  radial  components  respectively 
swirl  component  of  velocity. 


(4.  25) 

fd  z  (fd  2  "  id 

(4.  26) 

fdr  =  (fd2+fdl) 

2 

X 

2 

X 


V2sin0  , 
Vr  (l-cose) 


and  noting  that  the  mean  velocities  and  fluctuations  satisfy  these 
equations  independently, 

Sfdz  2 

(4.27)  =  |  sine  , 

z 

and 

3*  9 

(4.28)  |  (1-cose)  . 

r 

The  broadening  due  to  the  axial  fluctuations  is 

9  f 

ifz  -  ivi  STT-  =  !  aVzsine  ■ 


(4.  29) 
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which  for  this  experiment  was 


Af  =*  (10)  HZ. 
z 


Similarly,  the  broadening  due  to  the  radial  fluctuations  is 


(4.30) 

which  for  this  work  was 


Af  =  AV'  (1  -  cos  0) 
r  X  r 


Afy  «  (10)  HZ 

Thus  the  axial  turbulent  fluctuations  were  essentially  unaffected 
by  aperture  broadening,  whereas  the  radial  fluctuations  were  of  the 
same  order  of  magnitude  as  the  aperture  broadening. 

The  broadening  due  to  a  finite  size  measuring  volume  is 


(4.31) 


Af 


9 

•gA  AU 
3  V  z 


From  equation  (4 . 25) 
(4.32) 


af 


d  z 


3  V. 


3fdz  9Vz 
3Vz  3r 


sin  9  . 


Ar 


Here  A  r  is  just  the  effective  depth  of  focus  of  the  receiving  lens 
combination  of  and  Lr,  of  Figure  3. 1.  The  depth  of  focus  is  related 
to  the  wavelength  X  and  the  lens  number,  F,  by 


(4.33) 


Ar  =  4  FA 


where  the  lens  number  is  equal  to  the  lens  focal  length  divided  by 
the  aperture  size.  Therefore 


A  f  =  8  sin  9  F 


2  3Vz 
3  r 


(4.34) 
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Thus,  measuring  volume  broadening  becomes  strongest  in  the  region 
of  the  steepest  volocity  gradient  which,  for  this  experiment,  was  near 
the  wall. 

The  effect  of  finite  measuring  volume  size  can  be  minimized  by 
making  the  lens  number  of  the  receiving  lens  as  small  as  possible. 


Chapter  Y 
FLOW  SYSTEM 

L  Description  of  the  Flow  System 

The  flow  system  used  in  this  experiment  is  shown  in  Figure  5. 1. 
The  system  can  he  described  as  follows:  Water  from  the  city  water 
supply  entered  through  valve  and  passed  through  a  turbine  flow 
meter.  The  flow  rate  was  proportional  to  the  resulting  output  frequency 
which  was  monitored  on  a  frequency  counter. 

The  polymer  was  injected  into  an  elbow  near  the  mixing  tank 
inlet.  The  resulting  dilute  polymer  solution  passed  through  a  honey¬ 
comb  straightener  and  entered  a  two  inch  brass  pipe  via  a  rounded 
inlet  at  the  bottom  of  the  mixing  tank.  The  scattering  particles  were 
injected  just  beyond  this  inlet  by  a  piston  pump  which  was  driven  by 
air  pressure. 

The  pressure  drop  in  the  brass  pipe  was  measured  by  an  inclined 
U  -tube  manometer  attached  to  sidewall  pressure  taps  located  11.875 
feet  apart. 

The  fully  developed  flow  passed  through  the  test  section  (220 
diameters  from  the  pipe  inlet),  through  the  flow  control  valve,  V  , 
and  out  the  drain. 

A  single -pass  flow  system  was  chosen  rather  than  a  closed  loop 
system  to  avoid  mechanical  degradation  of  the  polymer  solution  and 
to  provide  a  freshly  mixed  solution. 
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2.  Polymer  Injection 

The  usual  method  of  preparing  a  polymer  solution  from  the  polymer 
powder  consists  of  preparing  a  master  batch  which  is  diluted  to  the 
desired  polymer  concentration.  Since  water  soluble  polymers  dissolve 
slowly  and  have  a  tendency  to  clump,  this  process  can  be  tedious  and 
time  consuming.  Furthermore,  since  the  solutions  are  usually  prepared 
well  in  advance  of  use,  degradation  due  to  aging  is  introduced. 

To  avoid  these  problems,  the  polymer  was  injected  as  a  powder 
into  the  water  flow.  This  injection  was  accomplished  by  pumping  a 
slurry  of  polymer  powder  in  methanol  into  the  water  flow.  Uniform 
concentration  of  the  polymer- methanol  solution  was  assured  by  con¬ 
tinuously  stirring  it  in  the  reservoir.  Since  methanol  is  a  non-solvent 
of  Separan  AP30  (the  polymer  used),  there  was  no  chance  for  clumping 
to  occur.  The  pump  used  was  a  variable  speed  tubing  pump  that  was 
regulated  to  compensate  for  changes  in  torque  or  system  pressure 
to  assure  a  constant  flow  rate , 

The  resulting  polymer-water  solution  was  quite  fresh  albeit  perhaps 
of  non-uniform  concentration  due  to  the  presence  of  macroscopic  lumps 
of  high  concentration  polymer.  It  was  hoped  that  this  non -homogeneous 
concentration  would  produce  a  variation  of  the  drag  reduction  phenomena 
usually  observed. 

The  polymer  used  was  a  polyacrylamide  (Separan  AP30).  This  type 
of  polymer  is  a  good  drag-reducing  agent,  and  is  also  less  susceptible 


to  shear  degradation  than  polyethylene  oxides,  which  are  also  good 
drag  reducing  agents. 

3.  Scattering  Particles 

Particle  contaminants  must  be  present  in  the  flow  in  order  to 
generate  a  Doppler-shifted  light  wave  through  light  scattering.  Yeh 
and  Cummins  (1964),  in  their  pioneering  work  on  the  LDA,  showed 
that  coherent  mixing  is  possible  using  light  scattered  from  uniform 
sized  polystyrene  particles.  Later,  Foreman  et  al  (1965)  showed 
that  optical  heterodyning  is  possible  using  scattering  particles  of 
irregular  size  and  shape. 

Natural  contaminants  cause  sufficient  scattering  for  heterodyning, 
when  using  high-power  lasers.  Since  a  5  mw  laser  was  used  in  this 
work,  it  was  decided  to  use  polystyrene  particles  of  uniform  size 
(average  diameter  0.  2  microns)  as  the  scattering  medium.  These 
particles  are  sufficiently  small  as  well  as  being  practically  neutrally 
buoyant,  so  that  they  do  not  alter  the  flow  (Huffaker,  1970).  This 
is  important  since  the  velocity  measured  by  the  LDA  is  the  scattering 
particle  velocity.  The  concentration  in  the  scattering  volume  was 
high  enough  to  prevent  a  significant  drop-out  in  the  heterodyned 
signal. 

The  cost  of  commercially  available  polystyrene  particles  in  the 
quantities  needed  is  prohibitively  high.  Therefore,  it  was  decided  to 
manufacture  them  from  raw  styrene  using  a  technique  developed  by 
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Woods,  Dodge,  and  Kriger  (1988).  The  particles  so  produced 
have  been  found  suitable  for  our  purposes. 

The  scattering  particles  were  injected  into  the  flow  just  below  the 
accumulator  tank  using  a  piston  (see  Figure  5.1).  The  piston  was 
driven  using  air  pressure.  A  velocity  transducer  measured  the  piston 
speed,  thus  enabling  calculation  of  the  scattering  concentration 
provided  the  water  flow  rate  and  the  particle  concentration  in  the 
cylinder  were  known  also. 

The  optimal  concentration  was  determined  by  varying  the  particle 

flow  from  the  cylinder  and  observing  the  effect  on  the  quality  of  the 

heterodyned  signal  as  seen  on  an  oscilloscope.  By  this  means,  the 

-3 

particle  concentration  was  chosen  to  be  2(10)  %by  weight.  This  is 
the  same  concentration  used  by  Chu  and  Schoenes  (1968)  in  their 
work. 


Chapter  VI 


SIGNAL  PROCESSING 

1.  Readout  Unit 

To  obtain  the  mean  velocity  and  turbulence  intensity  profiles  of  the 
axial  and  radial  components  in  a  pipe  flow,  one  must  be  able  to  extract 
the  mean  frequency  and  the  fluctuations  about  it  from  the  time  varying 
Doppler  frequency  output  of  the  P.  M.  Several  investigators  have  used 
spectrum  analyzers  to  do  this  (Goldstein  and  Kried,  1967;  Rudd,  1972). 
Unfortunately,  one  cannot  obtain  the  energy  spectra  or  Reynolds' 
stresses  using  this  signal  processing  technique. 

Others  (Foreman  et  al,  1965;  Fridman,  Kinnard,  and  Meister, 

1969)  have  developed  a  variety  of  signal  processing  devices  to  extract 
the  mean  and  fluctuating  frequencies  from  the  P.M.  output  signal.  These 
units  have  the  advantage  of  yielding  an  output  proportional  to  the  instan¬ 
taneous  velocity  fluctuations,  thereby  allowing  processing  to  obtain 
energy  spectra  and  Reynolds'  stress. 

The  unit  used  in  this  work  was  a  model  MFT  100  X  readout  unit 
made  by  Electronetics,  Inc.  It  has  the  capability  of  measuring  the 
mean  Doppler  frequency  as  well  as  the  time  varying  fluctuations  about 
the  mean.  This  unit  has  a  frequency  response  of  50  KHZ  to  10  MHZ 
which  was  found  sufficient  for  this  experiment. 

Referring  to  Figure  6. 1,  the  readout  unit  circuitry  can  be  described 
as  follows:  the  output  of  the  P.M.  tube  (the  heterodyned  signal)  is  fed 
into  the  input  jack  of  the  readout  unit.  From  here  it  goes  into  a  low 
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Readout  Unit  Circuitry 
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pass  filter  (10  MHZ  cut-off)  and  is  then  amplified  in  a  wideband  video 
amplifier.  The  amplifier  output  goes  to  a  front  panel  jack  so  that  the 
signal  may  be  viewed  on  an  oscilloscope.  Its  output  also  goes  to  a 
mixer-tuner  where  it  is  heterodyned  with  the  unit  oscillator  signal. 

The  resulting  signal  is  fed  into  a  frequency -to-voltage  converter  which 
yields  a  voltage  proportional  to  the  difference  frequency.  After  several 
stages  of  amplification,  the  signal  is  fed  into  a  front  panel  jack,  and 
also  to  a  front  panel  null  meter.  The  unit  oscillator  frequency  is 
adjusted  until  the  meter  indicates  a  null.  This  frequency  represents 
the  mean  velocity.  The  output  signal  is  then  proportional  to  only  the 
turbulent  part  of  the  velocity. 

2.  Signal  Processing  Set-Up 

The  instantaneous  velocity  of  a  turbulent  flow  is  considered  to  be 
the  sum  of  a  time -invariant  mean  velocity  and  a  velocity  fluctuating 
about  this  mean,  that  is, 

(6. 1)  V  (t)  =  V  +  V*  (t)  . 

Recalling  the  frequency-velocity  relation  developed  in  Chapter  2 

(6.2)  fd  =  £  V(t)  -  <ks-k0). 

we  see  that  V  and  V'(t)  each  satisfy  the  above  equation  independently. 
Similarly  the  components  of  V  and  V'(t)  individually  satisfy  the  component 
frequency  velocity  relations 
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(6.3) 
and 

(6.4) 


vz(t>  -  r 

o 


f  _  f 

C  ,  dl  d2 


sin  e 


fd  1  +  *d  2 


vr<‘>  ■  r  • 

Keeping  these  equations  in  mind,  let  us  consider  the  signal  processing 
set-up  shown  in  Figure  6.  2.  The  output  of  the  P.M.  s  feed  into  their 
respective  readout  units.  Each  unit  oscillator  can  be  adjusted  so  as  to 
produce  a  null  as  indicated  by  a  D.C.  voltmeter  attached  to  the  output 
of  the  readout.  The  frequencies  of  the  unit  oscillators  at  this  setting 
are  the  mean  Doppler  frequencies  and  therefore  yield  the  mean  axial 
velocity  through  equation  (6.3).  The  fluctuating  voltages  from  the  unit 
outputs  are  added  and  subtracted  from  the  input  signals  to  yield  outputs 
related  to  the  fluctuating  components  of  velocity  (see  equations  (6. 3)  and 
(6.4)  ).  The  outputs  of  this  unit  can  be  taped  for  later  wave  analysis, 
fed  into  RMS  meters  to  obtain  turbulence  intensities,  and  also  multiplied 
and  read  on  a  D.C.  voltmeter,  the  result  being  proportional  to  the 
Reynolds’  stress. 


3.  Signal  Processing 

As  mentioned  in  the  previous  section,  the  components  of  V  and  V'(t) 
(the  mean  and  fluctuating  components  of  velocity,  respectively)  individ¬ 
ually  satisfy  equations  (6.3)  and  (6.4).  Since  the  detector  produces  only 
the  absolute  values  of  the  frequencies,  it  is  helpful  to  introduce  the 
quintities 
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a  =  (1  -  cos0)/sin  6  , 

SD  =  sign  (Vz  -  aVr)  , 

Sg  =  sign(Vz  +  aVr)  , 

SR  =  W 

cz  =  csd/2  fQ  sine  . 

Then 


(RMS)  values  for  the  fluctuating  velocities  V'  and  V' ,  i.e. ,  [V'  ] 
i  z  r  z 

2  ^ 

2111(1  ]  ,  These  turbulence  intensities  are  related  to  RMSj  and 

RMS,.,  the  meter  readings  of  the  RMS  meters  shown  in  Figure  6,  2, 

In  order  to  develop  these  relationships,  consider  the  signal  processing 
set-up  of  Figure  6.  2. 


The  fluctuating  D,C.  voltage  outputs  from  readout  units  1  and  2, 
Vj  and  V'2,  were  amplified  by  amounts  and  a2,  respectively.  Both 


tof^* 


Meter 


Figure  6. 2 

Signal  Processing  Set-Up 
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signals  were  fed  into  the  analog  unit  which  yielded  their  difference  and 
sum  as  outputs, 

(6.9)  D  =  a1V,1  -  a2V’2, 

(6.10)  S  =  ajV'j  +  a2V'2. 

where  D  denotes  the  difference  of  the  input  voltages  and  S  their  sum. 
Now  aj  and  a.^  are  chosen  such  that, 

(6. 11)  ^2^1  =  a1^2  ’ 


where  Sj  and  S2  are  the  frequency -to-voltage  sensitivities  in  megahertz 
per  volt  of  readout  units  1  and  2,  respectively. 

Multiplying  equation  (6.  9)  by 

(6.12)  SjD  =  SjajV’j  -  S1a2V’2. 

Using  equation  (6.  11) 


(6.13)  SXD  =  a1(S1V\)  -  a^S^). 

The  fluctuating  Doppler-shifted  frequencies  f^ 


to  and  V'2by 

(6. 14) 

-  sivl  - 

(6. 15) 

fd2  ■  S2V2- 

and  f^2  are  related 


Substituting  equations  (6. 14)  and  (6. 15)  into  equation  (6. 13) 


SjD  = 
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(6. 16) 


f1  f* 
ldl  xd2 


SXD 


Substituting  equation  (6. 18)  into  equation  (6. 7)  for  =  +  1.0  , 

(6.17) 


C  S, 

Vz  =  S  • 


Similarly 

(6. 18) 


v;  - 


CzSl 
a  a. 


D  . 


Equations  (6- 17)  and  (6. 18)  define  the  relationship  between  the  fluctuating 
components  of  velocity  and  the  outputs  of  the  sum  and  difference  analog 
circuit. 

The  analog  outputs  D  and  S  were  further  amplified  by  amounts  b^ 
and  b2,  respectively,  before  being  fed  into  RMS  meters  whose  readings 
were  RMS.^  and  RMSg. 

(6.19)  RMSj  =  bj  [D^j  2  , 

(6.20)  RMS2  =  b2  [S'  2]  2  . 


Therefore,  from  equation  (6. 17) 

— 5T  a  CzSl 

(6.2«  [Vf]  -  ^7  BMS2  • 

Likewise 

(6-22>  Ivr2l5  *  RMS1  • 

wnich  relates  the  RMS  meter  readings  to  the  turbulence  intensities. 
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The  Reynolds'  stress  is  defined  asV^V^,.  Now,  the  output  of 
the  product  and  averaging  analog  circuit  (PAAC)  is  related  to  the 
Reynolds’  stress.  The  inputs  to  the  PAAC  were 


(6.23) 

bf 

V  -  srrr  - 

Z  1 

(6.  24) 

b9  a.a 

V  ■  cl  • 

z  1 

The  output  of  the 

PAAC  was 

(6.25) 

PAAC  =  b1Sb2D 

aa,  b0 

(6.26) 

PAAC  ■  rrrf  VzVr  • 

Cz  S1 

Therefore, 

,  Cz  Sl2 

(6.27) 

V'  V'  »  ~~ - —  PAAC  . 

z  r  *i2»iV 

which  relates  the  Reynolds'  stress  to  the  output  of  the  PAAC. 

b,  a 

The  outputs  of  the  sum  and  difference  analog  circuit  ^—5 

2  1 

and  -- yy—g—  V’  were  taped.  They  were  subsequently  each  fed  into 
CzSl  r 

the  wave  analyzer  which  yielded  their  power  spectral  densities 

(6.28)  Fz(w)  =  kz  e1WT  V^(t)V^(t +“ r)  dr  , 

(6.29)  Fr(w)  =  kr  iy”  VJ  (t)  V;  (t  7  r)  dr. 

which  can  be  transformed  into  one -dimensional  energy  spectra  in  wave 
number  space  by  using  the  transformation 
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(6.30) 


kj  =  2ttw/Vz  , 


where  k,  is  the  wave  number.  The  wave  analyzer  was  calibrated  by 
lumping  all  of  the  calibration  factors  such  as  b^  a^,  b^,  etc.  into 
one  parameter  k.  This  parameter  was  determined  by 


(6.31) 


(6.32) 


T’z“  *  2  Jo  iy^w)dw’ 

vT2 

lr  -  _ ± _ 

z  r  oo  * 

2  4  Fz(w)  dw 


Likewise 


(6.33) 


k  = 
r 


2  4  dw 


where  V7  and  V7  were  known  from  previous  measurements  and  the 
z  r 

integrals  were  obtained  as  outputs  on  the  wave  analyzer. 


Chapter  VII 

MEASUREMENTS  AND  OBSERVATIONS 

1.  Discussion  of  the  Turbulence  Measuring  Technique 

The  stated  goals  of  this  work  were  to  measure  the  drag  reduction, 
mean  flow  characteristics,  and  turbulence  properties  in  the  axial 
and  radial  directions  for  drag  reducing  pipe  flows.  To  obtain  such 
data  with  the  beam  configuration  used  m  this  work,  one  must  generally 
process  the  signals  from  both  optical  units  simultaneously,  i.e. ,  the 
readout  unit  outputs  must  be  fed  into  the  analog  circuit  (Figure  6.  2) 
simultaneously.  It  is  not  possible  to  generate  meaningful  data  by 
recording  the  output  of  one  readout  unit  on  tape,  at  a  later  time  record 
the  output  of  the  other  readout  unit  on  tape,  and  to  then  read  both  tapes 
simultaneously  into  the  analog  circuit  of  Figure  6.2,  for  phase  rela¬ 
tionships  would  thereby  be  lost.  Concern  for  this  was  occasioned  be¬ 
cause  the  Electr one tics,  Inc.  readout  units  used  in  this  work  have  never 
been  in  working  order  at  the  same  time.  Thus  it  was  never  possible 
to  process  the  data  from  channels  1  and  2  simultaneously.  As  a  conse¬ 
quence  it  was  not  possible  to  measure  the  radial  intensities  or  V' V' 
(unaveraged).  What  was  done  instead  was  to  use  the  two-dimensional 
optical  geometry  to  measure  the  Reynolds'  stress  and  the  mean  velocity 
and  to  use  a  one -dimensional  geometry  (Chung  and  Graebel,  (1972) )  to 
measure  V^  and  its  spectra.  This  kept  the  loss  of  measurable  data  to 
a  minimum. 
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All  measurements  were  made  at  a  Reynolds’  number  of  50,000 
based  on  the  average  velocity  in  the  pipe.  Data  for  water  was  compared 
to  that  obtained  for  50,  100,  and  200  PPM  solutions  of  Separan  AP30 
and  water. 

As  mentioned  before,  the  polymer  injection  technique  used  in  this 
wurk  is  unique.  It  offers  many  advantages  over  the  master  batch 
technique  commonly  used  in  experimental  work.  As  such,  it  is  a 
promising  technique  for  large-scale,  accurate  concentration  polymer 
injection  in  laboratory  or  non-laboratory  applications.  It  was  expected 
that  this  injection  technique  would  show  significantly  different  drag 
reduction,  mean  flow,  and  turbulence  properties  than  observed 
previously.  That  such  was  the  case  is  shown  in  the  following  data. 

2.  Discussion  of  Results 
2. 1  Drag  Reduction 

Drag  reduction  was  measured  as  a  pressure  loss  by  means  of  two 
side  wall  pressure  taps  located  11.  875  feet  apart.  The  pressure 
difference  between  them  was  measured  on  a  U-tube  manometer.  In 
measuring  the  pressure  loss  of  the  polymer  solutions,  it  was  found 
that  maximum  drag  reduction  was  established  two  minutes  after  the 
polymer  injection  pump  was  started.  Therefore,  ali  flow  and  turbulence 
measurements  were  made  only  after  the  pump  had  been  operating  for 
two  minutes.  Drag  reduction  effects  were  noted  in  the  pipe  for  up  to 
eight,  minutes  after  the  pump  was  stopped.  These  time  lags  were  due 


to  he  presence  of  the  mixing  tank  upstream  of  the  pressure  probes. 


Tb«  percentage  drag  reduction  was  calculated  by: 
Percentage  Drag  Reduction  = 


AP  -  AP 

S  ^  x  100 


AP. 


where:  AP  :  pressure  loss  due  to  friction,  solvent  only, 

s 

APp:  pressure  loss  due  to  friction,  polymer  solution. 


The  percentage  drag  reduction,  friction  velocity  V  ,  the  ratio  of 
maximum  to  average  velocity,  and  friction  factor,  f,  for  various 
polymer  concentrations  are  presented  in  Table  7.  1.  The  friction 
factors  of  this  work  are  compared  to  those  of  Chung  and  Graebel  (1972), 
for  the  same  polymer,  polymer  concentration,  and  Reynolds’  number. 
Whereas  the  friction  factor  for  water  was  found  to  be  slightly  lower  in 
this  work,  the  friction  factors  for  the  polymer  solutions  are  higher 
than  found  by  Chung  and  Graebel  (1972).  This  indicates  that  less  drag 
reduction  has  occurred  in  the  present  study.  One  reason  for  this  might 
be  that  the  polymer  powder  was  not  completely  dissolved  and  therefore 
the  actual  polymer  concentration  was  less.  Since  drag  reduction 
increases  with  increased  polymer  concentration,  up  to  a  limit  (Virk. 
Mickley,  and  Smith,  1970), this  would  explain  the  lower  drag  reduction 
observed.  Solid,  undissolved  polymer  particles  were,  in  fact, 
observed  passing  through  the  test  section. 


The  ratio  of  maximum  to  average  velocities  were  higher  for  all  of 
the  polymer  concentrations  than  for  water.  This  contradicts  the  findings 
of  Chung  and  Graebel  (1972)  and  Seyer  and  Metzner  (1969).  The  reason 
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for  i.his  is  that  the  presence  of  the  solid,  undissolved  pol/mer  particles 
sign  ficantly  effected  the  mean  flow  and  turbulence  properties  of  this 
pipe  flow.  Subsequent  data  will  bear  this  out. 

2.2  Mean  Velocity  Profiles 

The  axial  mean  velocity  profiles  for  water,  50,  100,  and  200  PPM 
solutions  of  Separan  AP30  flowing  at  a  Reynolds'  number  of  50,000  are 
plotted  in  Figure  7. 1.  The  mean  velocities  have  been  normalized  by  the 
maximum  velocity,  i.e.,  the  centerline  velocity.  The  at  scissa  is  the 
normalized  distance  from  the  wall,  r  =  1  -  —  ,  where  r  is  the  distance 
from  the  flow  axis  and  a  the  radius  of  the  pipe.  The  confidence  band  of 
this  data  is  5  per  cent.  The  mean  velocity  profile  of  Sandborn  (1955),  who 
measured  air  flow  in  a  pipe  at  a  Reynolds'  number  of  50,000  is  also  shown 
in  Figure  7. 1,  His  results  are  seen  to  be  in  general  agreement  with 
those  of  the  present  work.  It  is  seen  from  Figure  7. 1  that  the  velocity 
profiles  tend  to  become  less  blunt  as  the  polymer  concentration  increases. 
This  is  in  direct  contradiction  with  the  work  of  Chung  and  Graebel  (1972) 
and  Rudd  (1972)  who  show  lunter  velocity  profiles  for  polymer  solutions 
than  for  water.  Both  studies  were  done  using  the  same  polymer  and 
similar  polymer  concentrations  at  the  same  Reynolds'  number.  This 
apparent  anomaly  is  just  a  manifestation  of  the  particulate  nature  of 
the  flow.  Daily  and  Hardison  (1964)  found  that  the  velocity  profile  for 
a  2')%  (by  volume)  concentration  of  spherical  plastic  particles,  0.004 
to  0.006  inches  in  diameter,  in  a  100,000  Reynolds'  number  flow  was 
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mu'*h  less  blunt  than  for  a  purely  water  flow.  Furthermore,  they 
fou.d  that  the  ratio  of  maximum  to  average  velocity  for  the  particulate 
flow  was  higher  than  for  the  water  flow.  This  correlates  with  what 
was  observed  in  the  present  work.  It  is  reasonable  to  conclude  that, 
for  the  core  region,  the  mean  flow  characteristics  of  the  polymer 
solution  are  dominated  by  particulate  effects. 

2. 3  Turbulence  Intensities 

The  axial  turbulence  intensity  profiles  are  shown  in  Figure  7.  2. 

The  intensities  have  been  normalized  with  respect  to  the  maximum  mean 
velocities.  As  mentioned  before,  the  intensity  profiles  were  measured 
using  a  one-dimensional  set-up  similar  to  the  one  used  in  the  work  of 
Chung  and  Graebel  (1972).  The  transformation  from  the  two-dimensional 
to  one -dimensional  set-up  was  accomplished  by  simply  replacing  B.  S.  2 
with  a  mirror.  It  was  hoped  that  an  included  beam  angle  (the  angle 
between  the  incident  and  reference  beams)  of  20  degrees  could  be  used, 
allowing  for  a  greater  instrument  sensitivity.  This  was  not  possible, 
however,  since  the  heterodyned  signal  at  this  angle  was  too  weak  and 
broadened  to  yield  meaningful  data.  The  weakness  of  the  signal  was 
due  to  the  low  forward  scattering  intensity  at  this  relatively  large 
included  beam  angle.  (Forward  scattering  intensity  decreases  rapidly 
with  increasing  included  beam  angle.)  The  signal  broadening  was  due 
to  he  large  measuring  volume,  also  a  consequence  of  tht  large  angle 
(Jackson  and  Paul,  1970).  Experience  has  shown  that  th<>  optimal 


O  Water 

•  50  PPM 

6  100  PPM 

□  200  PPM 

7  Sandbom  (1955) 
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included  beam  angle  is  approximately  12  degrees  for  this  set-up. 

Thus,  when  the  beam  angle  was  decreased  to  12.  4  degrees,  a  clear 
heterodyned  signal  was  observed. 

Referring  to  Figure  7.2,  it  is  seen  that  the  axial  turbulence 
intensities  are  everywhere  higher  for  the  polymer  solutions  than  for 
water.  However,  others  have  found  that  the  axial  intensities  are  lower 
for  polymer  solutions  than  for  water  in  the  range  0.4  <  r  <  1.0  and 
higher  in  the  range  0  <  r  <0.4  (Chung  and  Graebel,  1972)  and  (Rudd, 
1972).  That  the  intensities  are  higher  for  the  polymer  solutions  in  the 
region  0.4  <  r  <  1. 0  in  this  work  is  a  consequence  of  the  particulate 
nature  of  the  flow  in  the  core  region.  The  particles  in  the  flow  generate 
eddies  which  are  manifested  as  higher  turbulence  intensities.  This 
was  found  experimentally  by  Daily  and  Hardison  (1964)  who  noted  higher 
axial  intensities  for  a  particulate  flow  than  for  water,  in  a  pipe  flow. 

The  intensities  of  Figure  7.  2  show  a  marked  decrease,  all  across  the 
profile,  with  increasing  polymer  concentration.  Thus,  the  trend  of 
the  polymer  intensities  is  toward  the  usually  observed  drag  reducing 
polymer  effect.  This  would  seem  to  indicate  that  the  effect  of  the 
increase  of  polymer  concentration  is  stronger  than  the  effect  of  the 
increase  in  particle  concentration.  It  would  seem  that  at  some  higher 
polymer  concentration  than  measured,  the  polymeric  action  would 
dominate  the  particulate  effect  to  the  extent  of  showing  the  usually 
observed  turbulence  properties  of  drag  reducing  flows.  Subsequent 
Reyr olds'  stress  and  energy  spectra  data  supports  this  hypothesis. 
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2.4  Reynolds*  Stress 


The  normalized  Reynolds'  stress,  V'  V'  /V  ,  is  shown  in 
Figure  7.3  plotted  against  r.  The  Reynolds'  stress  can  be  thought  of 
as  a  measure  of  the  rate  of  transport  of  z-  momentum  through  a  surface 
normal  to  the  r-axis  by  the  turbulent  fluctuations.  It  may  also  be 
thought  of  as  an  additional  stress  due  to  the  presence  of  turbulent 
velocity  fluctuations.  In  light  of  this,  it  would  be  expected,  that  for  a 
particulate  flow  (at  least  in  the  core  region)  the  Reynolds'  stress 
would  be  higher  than  for  the  solvent  flow.  That  the  Reynolds'  stress 
for  the  polymer  solutions  are  higher  than  for  water  reiterates  the 
particulate  nature  of  the  polymer  flow. 

The  turbulent  energy  equation  for  fully  developed,  incompressible, 
axially  symmetric  flow  is: 
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where  the  term  -g— 


represents  a  net  rate  of  exchange  of  energy  between  the  mean  flow  and 
the  turbulent  field  through  the  working  of  the  Reynolds'  stress  against 
the  mean  velocity  gradient.  In  other  words,  it  represents  the  production 
of  turbulent  energy.  It  has  been  postulated  (Walsh,  1967)  that  the  small 
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eddies  near  the  wall  which  grow  and  form  the  large-scale  turbulence, 
gain  energy  through  interaction  with  the  mean  velocity  profile  through 
its'  Reynolds’  stress  and  lose  energy  by  viscous  dissipation  and  local 
storage  in  the  polymer  molecules.  If  the  eddies  gain  more  energy  than 
they  lose,  they  tend  to  grow.  If  they  lose  more  energy  than  they  gain, 
they  will  be  annhilated  by  viscous  dissipation.  The  effect  of  the  polymer 
is  to  inhibit  growth  of  these  eddies  and  to  facilitate  their  destruction 
through  the  action  of  viscosity.  Thus,  there  will  be  fewer  small  scale 
eddies  per  unit  volume  per  unit  time  and  as  a  consequence,  less  large 
scale  turbulence.  This,  Walsh  feels,  is  the  mechanism  for  drag 
reduction.  In  light  of  this  theory,  it  would  be  expected  that  the  pro¬ 
duction  term  should  be  less  for  the  polymer  solution  than  for  the  solvent, 

at  least  near  the  wall.  Refering  to  Figure  7. 1  and  7. 3,  it  is  seen  that 
dV 

V’ Vr  — s —  is  less  for  the  200  PPM  polymer  solution  than  for  water  in 
z  t  u  r 

the  region  near  the  wall.  This  would  seem  to  indicate  that,  for  the  200 
PPM  solution,  the  polymeric  action  is  dominating  the  particulate  effect 
in  the  turbulent  properties. 

2.  5  Energy  Spectra 

One-dimensional  energy  spectra  are  plotted  versus  one -dimensional 
wave  number  space  in  Figures  7.4,  7.5,  and  7.6  for  various  polymer 
concentrations  at  three  radial  locations;  r  =  1.00,  0.500,  and  0.020, 
respectively.  These  spectra  were  obtained  b„  tape  recording  the  axial 
velocity  fluctuations  (the  bandwidth  of  the  tape  recorder  being  D.  C.  to 
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2500  Hz)  and  later  reading  this  taped  data  into  a  wave  analyzer.  The 
data  was  analyzed  at  a  13  Hz  sec  sweep  rate  using  a  filter  whose 
effective  bandwidth  was  13  Hz.  The  sample  time  was  180  seconds.  Thus, 
the  expected  error  is  about  2%. 

There  are  no  valid  conclusive  energy  spectra  of  polymer  flows. 

Of  the  few  investigators  who  have  attempted  energy  spectra,  the  validity 
oi  -.he  work  of  Virk  et  al  (1967),  Wells,  Harkness  and  Meyer  (1968),  and 
Johnson  and  Barchi  (1968)  must  be  seriously  questioned  due  to  the  basic 
measuring  techniques  that  they  employed.  (This  was  discussed  in  detail 
in  Chapter  1).  The  energy  spectra  of  Chw^g  and  Graebel  (1972),  who 
used  a  LDA  hi  their  work,  are  inconclusive  as  they  show  no  general 
trend.  Others  who  have  studied  drag  reducing  flows  with  laser  anemo- 
,0.9 ter s  have  not  reported  energy  spectra  because  diey  processed  their 
P.M.  outputs  with  spectrum  analyzers.  Thus,  the  time  domain  is 
annhilated  and  there  can  be  no  transformation  to  the  frequency  or  wave 
number  domain  required  to  produce  energy  spectra. 

Referring  to  Figures  7.  4,  7.  5,  and  7.  6,  it  is  seen  that  the  polymer 
spectra  ar-.  generally  higher  in  the  low  wave  number  (energy  containing) 
range  at  all  radial  positions  than  the  water  spectra.  A.s  it  would  be 
expected  that  a  particulate  flow  would  be  r.  energetic  than  a  solvent 
flow,  this  is  another  reflection  of  the  pi-.r.ieulate  nature  of  the  polymer 
Hows,  however,  the  energy  spectra  decrease  with  increasing  polymer 
concent  rat  on  in  this  wave  number  range.  Since  it  is  expected  that 
polymeric  action  lowers  the  energy  spectra  profile,  this  again  shows 


the  dominance  of  the  polymeric  effect  over  the  particulate  effect.  The 
spectra  do  not  generally  exhibit  the  -5/3  slope  predicted  by  Kolomogoroff's 
hypothesis.  This  is  not  surprising  since  the  Reynolds’  number  for  these 
flows  was  only  50,000  and  therefore  probably  not  high  enough  to  satisfy 
the  conditions  of  this  hypothesis  (Hinze,  1959).  Of  special  interest  are 
the  spectra  at  r=0.020  (Figure  7.  6)  which  is  near  the  wall.  It  is  seen  that 
spectra  for  the  100  and  200  PPM  concentrations  are  less  in  the  high  wave 
number  (the  dissipative)  range  than  the  water  spectra.  This  confirms 
Walshs'  theory  that  there  are  fewer  small  eddies  per  unit  volume  per 
unit  time  in  polymer  flows  than  in  water  flows.  Furthermore,  the  spectra 
of  the  200  PPM  solution  is  everywhere  less  than  the  water  curve,  showing 
the  complete  dominance  of  the  polymeric  action  in  the  wall  region  at 
this  concentration. 


3.  Summary  and  Conclusions 

The  drag  reduction,  mean  flow,  and  turbulence  properties  of  drag 
reducing  flows  were  measured  using  a  two-dimensional  laser  anemometer 
technique. 


The  optical  bench  used  was  designed  to  offer  experimental 
flexibility  ard  vase  of  alignment  as  well  as  to  satisfy  the  spatial  require¬ 
ment''  of  photomixing  and  to  minimize  spectral  broadening  inherent  in 
the  optical  set-up.  The  signal-to-noise  ratio  of  the  optical  set-up  and 
th  ?  photomultiplier  tube  was  studied  and  the  system  was  operated  in 


oider  to  maximize  this  ratio. 
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A  flow-through  system  was  chosen  for  this  work  to  minimize  the 
effects  of  polymer  degradation.  The  polymer  injection  system  was 
unique  to  laboratory  work.  The  polymer  powder  was  injected  directly 
into  the  water  flow  as  a  slurry  of  polymer  and  methanol.  This  technique 
was  used  to  prevent  both  polymer  aging  and  the  tedium  associated  with 
the  commonly  used  master-batch  technique. 

The  key  to  the  signal  processing  set-up  is  the  readout  unit  which 
takes  the  P.M.  signal  and  processes  it  to  yield  the  mean  velocity  and 
a  D.  C.  ijHignal  directly  proportional  to  the  fluctuating  component  of 
velocity.  Two  of  these  units  operating  simultaneously  with  a  two  - 
dimensional  LDA  allow  signal  processing  to  be  accomplished  with  the 
same  equipment  used  to  process  signals  from  X-wires  in  hot-wire 
anemometry.  Unfortunately,  the  readout  units  have  not  reached  the 
stage  of  development  where  they  are  completely  reliable.  In  fact,  the 
two  units  used  in  this  work  were  never  both  in  operating  order  simul¬ 
taneously.  Thus,  it  was  not  possible  to  make  all  of  the  measurements 
potentially  possible  with  a  two-dimentional  LDA. 

The  data  generally  contradict  the  work  of  previous  investigators. 
The  percentage  drag  reduction  for  50,  100,  and  200  PPM  solutions  of 
Separan  AP30  and  wate-  is  lower  than  that  found  previously.  This  is 
the  first  clue  that  the  polymer  powder  injected  did  not  dissolve  com¬ 
pletely.  The  mean  velocity  profiles  were  less  blunt  and  the  ratio  of 
maximum  to  average  velocities  were  higher  for  the  polymer  solutions 
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than  for  water.  This  disagrees  with  previous  work  and  is  an  indication 
of  he  particulate  (undissolved  polymer  powder)  nature  of  the  flow.  The 
turbulence  intensities  and  Reynolds'  stresses  are  higher  for  the  polymer 
solutions  than  for  water,  at  least  in  the  core  region.  This  also  is  not 
what  would  be  expected  for  the  usual  drag  reducing  flow  but  rather  for 
a  particulate  flow.  The  higher  energy  spectra  observed  for  the  polymer 
solutions  than  for  the  water  in  the  energy  (low)  wave  number  range  is 
another  indication  of  the  essentially  particulate  nature  of  the  polymer 
flows. 

Thus  the  data  shows  that,  for  low  polymer  concentrations,  using 
this  polymer  injection  technique,  the  flow  is  dominated  by  particulate 
effects.  As  polymer  concentration  increases,  however,  polymeric 
effects  are  seen  to  dominate  the  particulate  effects  to  the  extent  that  at 
the  200  PPM  concentrations  the  usual  drag  reduction  phenomena  begins 
to  appear  in  the  turbulence  properties.  This  is  especially  noticable  in 
the  wall  region  where  polymeric  action  is  the  strongest. 
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Appendix  A 


REFRACTIVE  EFFECT  OF  TEST  SECTION  AND  WATER 
ON  BEAM  GEOMETRY  OF  A 
TWO-COMPOSENT  LASER-D0PPLER  ANEMOMETER 


The  velocity  to  frequency  relationship  was  shown  in  Chapter  II 


to  be 
(A.  1) 
and 
(A.  2) 


Vz  = 


V£d2  "  fdl^ 

2  Bind 


y  =  Xe  ^fd2  +  fdP 

r  2(1 -cosd) 


where  X^  is  the  wavelength  of  the  incident  light  in  the  fluid  flow 
and  6  is  the  angle  of  incidence  in  the  fluid  flow.  Due  to  the  refractive 
properties  of  the  test  section  and  the  water  flow,  this  will  be 
different  than  its  value  in  air.  Likewise  the  angle  of  incidence,  6, 
will  be  different  than  the  angular  position,  cp,  (see  Figure  A.  1)  of  the 
reference  beams  before  they  hit  the  test  section. 

The  relationship  between  0  and  $ can.be  found  by  considering 
Figure  A.  1  and  Snell's  law: 

sin  Q 


(A.  3) 


Q 


sin  <fr 

n  , 

<P 


where  n  is  the  index  of  refraction. 
Thus 


(A.  4) 


-1  nfi 

6  =  sin  [  — —  sin  $] 
0 
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If  we  note  that 


X  :  wavelength  of  light  in  the  flow 

U 

X  :  wavelength  of  light  in  a  vacuum 


equation  (A.  1)  can  be  written  as 

l0X  (fd  2 


(A.  5) 


Vz  = 


2  sind 


But  in  light  of  equation  (A.  3)  this  becomes 

n0X  (fd  2  "  W 


(A.  6) 


Vz  * 


2  sin  0 


or  since  n  =  1 

(A.  7)  Vz  =  X(fd2-fdl)/2sin0  . 

Thus  the  form  of  equation  (A.  1)  remains  the  same.  The  form 
of  equation  (A. 2)  would  change,  however,  if  we  were  to  write  it  in 
terms  of  the  air  variables  X  and  0. 

Let  us  next  consider  the  following  problem:  given  that  the  beams 
are  moved  an  external  distance  dx  (see  Figure  A.l),  what  is  the  dis¬ 
placement  of  their  intersection  dr  ? 

Referring  to  Figure  A.  1  , 


(A.  8) 


y  =  r  tan  Q  +  t  tan  0  +  x  tan  0  , 


(A.  9) 


0  =  <ir  tan  8  +  dx  tan  0  . 


For  small  angles 


(A.  10)  tan  6  «  sin  9 

and 

(A.  11)  tan  <b  »  s in  0 


Therefore 
(A.  12) 


But  from  equation  (A.  3) 
(A.  13)  dr 


-n  . 

= 

n9 


which  for  our  set-up  becomes 


Appendix  B 


MEASUREMENT  TECHNIQUE  FOR  MEASURING 
THE  SWIRL  COMPONENT  OF  VELOCITY 

1.  Frequency-to-Velocity  Relationship 

The  optical  bench  used  in  the  present  work  has  the  flexibility  to 
alio  measurement  of  the  third  (swirl)  component  of  velocity.  To 
measure  the  swirl  component,  an  incident  beam,  k^,  (see  Figure 
B.  1)  is  aligned  perpendicular  to  the  flow  axis  at  an  angle  0  with  respect 
to  the  X  axis.  Incident  beam  k  of  equal  intensity,  (not  shown  in 
Figure  B.  1)  is  symmetrical  to  k^  about  the  X  axis.  The  beams  inter¬ 
sect  at  a  distance  6  from  the  flow  axis.  Thus,  beams  and  kQ2 
form  a  plane  perpendicular  to  the  flow  axis  and  are  symmetrical  about 
the  X  axis.  Recalling  the  frequency -to-velocity  relationship  developed 
in  Chapter  2 

<B-‘>  fd  -  r  X  -  fes-y 

where  for  this  beam  geometry 


(B.  2) 

V  -  V  e  + 
—  r  —  r 

V*’ 

(B.3) 

k„,  =  cos0e  + 

-01  —  r 

sin9£$  > 

(B.4) 

*02  =  c°s9Sr  ' 

sinee0  , 

and 

(B.  5) 

k  =  k  1  = 

-s  -si 

—  s2  ' 
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Thi?  last  equation  implies  that  we  are  sensing  (with  a  P.M.)  the 
scattering  due  to  bothk^  and  along  the  same  direction,  kg  . 
Continuing, 

f 

(Be  8)  fdl  =  £  [  -  Vr  cos  Q  -  Vg  sin  0  +  V  •  kg  ] 

and 


(B.  7) 


d2 


o 

C 


[-V  cos#  +  V\sine  +  V 
1  r  o  — 


k  1 

—  s  1 


As  seen  in  Chapter  IV,  the  P.M.  tube  senses  the  difference  frequency 
of  two  beams  impinging  on  its  face  which,  in  this  case,  is  f  d2  -  f'd  ^  . 


(B.  8) 
(B.9) 


fd  2  “  fdl  =  2  C  V0Sin  6  * 


V  =  C  2  ~  W 
6  f  2  sin  Q 


Ijet  us  consider  this  last  equation.  It  tells  us  that  the  frequency- to- 
velocity  relation  is  not  a  function  of  scattered  beam  position  and  there¬ 
fore  is  not  a  function  of  P.M.  position.  As  a  consequence,  there  is  no 
aperture  broadening.  This  beam  geometry  and  its  advantages  were 
first  noted  by  Mezumder  and  Wankum  (1970)  who  call  it  the  differential 
Doppler  anemometer. 


2.  Measurement  of  Turbulent  Fluctuations  About  a  Zero  Mean 

In  the  pipe  flow  being  studied,  the  only  non-zero  mean  velocity 
is  in  the  axial  direction.  Thus  the  radial  and  swirl  velocities  have 
only  fluctuating  components 
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(B.  10) 

vr(t)  =  v;  (t) 

and 

(B.  11) 

V‘>  *  v'e  «  ' 

However,  t}ie  bandwidth  of  the  readout  unit  used  in  this  work  is 
50  KHZ  -  10  MHZ.  It  will  not  measure  mean  velocities  or  fluctuations 
about  the  mean  for  velocities  corresponding  to  frequencies  under 
50  KHZ.  Therefore,  the  readout  unit  cannot  measure  V*  using  the 

C7 

beam  geometry  outlined  above  unless  one  creates  an  artifical  mean 
velocity.  This  can  be  accomplished  by  shifting  the  frequency  of  one  of 
the  incident  beams  using  a  rotating  radial  diffraction  grating  (Denison 
and  Stevenson,  1970).  Using  this  device  one  can  shift  the  incident 
frequency  by  AfQ,  where 

(B.  12)  M  =  w/K  . 


where  w  is  the  angular  velocity  of  the  radial  diffraction  grating  and  K 
is  a  grating  constant.  Frequency  shifts  of  up  to  500  KHZ  can  be  easily 
effected  using  thi6  technique. 

Let  us  assume  that  the  frequency  of  the  incident  beam  of  Figure 
A.  2  is  shifted  by  Af  j.  The  scattering  frequency  will  be 

(fo  +  Afl) 

teas)  fsl  =  VAfi+ — c — X  *  ~  — oi}  • 


(B.  14) 
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Note-  that  k  =  k  ,  =  k  0  since  a  differential  Doppler  configuration 

—  S  — SI  S  Li 


is  being  used.  The  P.M.  tube  sees  the  difference  frequency 
03.  IS) 


fsl  "  fs2 


_  Afl  +  C  -  ^02  "  -01'  +  C  ~ 


If  f  j  -  f  is  to  take  on  a  simple  form, 


(B.  16) 


Ss  ’  ioi  • 


Therefore 

(B*  I?)  fD  =  fsi-fs2  =  Afl  +  CT  -  ^—02  "  -01^  * 

which  becomes,  refering  to  Figure  A.  2 

2f 

(B.  18)  fD  =  Afx  +-^V,Qsin9  „ 

Thus,  Afj  represents  a  mean  velocity  and  the  second  term  the  frequency 
fluctuating  about  Afj.  Note,  however,  in  simplifying  the  frequency- 
velocity  relationship  the  advantage  of  the  differential  Doppler  has 
been  lost  since  the  position  of  k  g  has  been  defined. 

3.  Refractive  Effect  of  Test  Section  and  Water  on  Beam  Geometry 
In  order  to  measure  a  velocity  profile  across  the  test  section, 
it  is  necessary  to  know  the  relationship  between  the  incremental  move¬ 
ment  of  the  optical  bench  and  the  resulting  incremental  motion  of  the 
beam  intersection  (measuring  volume)  within  the  test  section.  This 
relationship  is  not  simple  due  to  the  refraction  of  the  incident  beam 
by  the  test  section. 


Assume  the  initial  distance  of  the  incident  beam  from  the 


reference  beam  at  the  front  surface  of  the  test  section  is  known  to  be 
y^.  If  the  optical  bench  is  moved  a  distance  Ax  in  the  negative  X 
direction,  the  incident  beam  intersects  the  test  section  a  distance 
y  =  y.  +  Axtan^from  the  reference  beam.  Using  Snell's  law  (see 
Figure  R.  1) 

a 

(B.19)  ip  =  arcsin(— ^ -  sin$  . 

<t> 

Noting  that  the  distance  (a  +  t)  from  the  center  of  the  test  section  to 
its  front  edge  can  be  written  as 

(B.20)  acose  +  bcos tp  =  a  +  t  , 

also 

(B.  21)  a  sine  +  bsin  ip  =  y  . 

These  equations  can  be  combined  to  yield 

(B.22)  sin(0  -  ip)  =  simp  -  2.  cos  \p  . 

ft  ft 

Again  using  Snell's  law 

nn 

(B.  23)  sine  =  ~  sin  {e  -  \p) 

V 

or 

(B.  24)  sine  =  ~  ~ -  sin  <p-  -~cos(sin  )] 

ip  (p  <p 

Now 


:b.2 
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(B.  26) 

and 
(B.  27) 


a  sin  9 
sin  9  +  e 


6  =  a  cos  9  -  S.  cos  (6  +e)  . 


Equations  (B.22),  (B.24),  (B.26),  and  (B.27)  can  be  solved  successively 
to  obtain  6,  the  incremental  movement  of  the  beam  intersection  relative 


